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Abstract 

The Dualized Standard Model has scored a number of successes in 
explaining the fermion mass hierarchy and mixing pattern. This note 
contains updates to those results including (a) an improved treat- 
ment of neutrino oscillation free from previous assumptions on neu- 
trino masses, and hence admitting now the preferred LMA solution to 
solar neutrinos, (b) an understanding of the limitation of the 1-loop 
calculation so far performed, thus explaining the two previous discrep- 
ancies with data, and (c) an analytic derivation and confirmation of 
the numerical results previously obtained. 



1 Introduction 



The Dualized Standard Model (DSM) suggested a few years ago [H 121 El has 
scored some, to us, notable successes in explaining the intricate and otherwise 
mysterious mass and mixing patterns of quarks and leptons. With only 3 
adjustable parameters, the model is able to reproduce already at 1-loop level 
the following mass and mixing parameters all to within present experimental 
errors |H EJ: all 9 CKM matrix elements \V rs \ for quarks the 2 MNS 
lepton mixing matrix elements | and \ U e s\ bounded by neutrino oscillation 
experiments with respectively atmospheric [7] and reactor neutrinos [S], plus 
the 3 mass ratios m c /mt,m s /mb and m M /m r [6 , together accounting for 8 
independent parameters of the Standard Model. However, the treatment so 
far published, culminating in (5] to which we refer the reader for details, 
contains some shortcomings which were not apparent at the beginning but 
have since, with time, been recognized and remedied. In particular, (a) 
some assumptions were made in v oscillations which are later found to be 
unnecessary, (b) the limitations of the 1-loop approximation made which 
were not recognized before are now seen to explain why the other masses 
and mixing angles: m u , m e , U e 2, were not well reproduced, while (c) the 
results obtained only numerically before are now seen to follow also from 
simple analytic considerations. The purpose of this note is to remove these 
old shortcomings, resulting in both a tighter and a more transparent scheme, 
while affording at the same time an analytic check on the previous numerical 
results. 

We begin with a very brief summary of the basic features of the DSM 
scheme necessary for our presentation. For explanations and details, the 
reader is referred to the literature cited. The DSM scheme is based on a 
nonabelian generalization of electric-magnetic duality which offers an ex- 
planation for the existence of exactly 3 fermion generations and a suggestion 
of how the generation symmetry is broken. This leads in turn to the con- 
struction of a Higgs potential and Yukawa coupling which give the tree-level 
fermion mass matrix a factorized form pQ and allows its radiative correc- 
tions to be calculated. Further, the fermion mass matrix is found to retain 
a factorized form even after 1-loop radiative corrections: 



where only the normalization ttit depends on the fermion species, i. e. U or 
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D type quarks, or charged leptons L or neutrinos N. The vector (x,y,z), 
however, now rotates in generation space with changing scales /i, and satisfies 
the RG equation: 



(2) 



with 

Si = «(^) + «^) cyclic> (3) 
x A + y A x z + z z 

and p 2 the Yukawa coupling strength |lj . In (J2J), one sees that the vector 
(x,y,z), taken to be normalized from now on, is stationary in orientation 
with respect to change in \i at (x,y,z) = (1,0,0) and ^(1,1,1), which 
means that these are to be interpreted as fixed points of rotation. Further, 
the sign of the derivative is such that as fi changes, the vector r = (x, y, z) 
traces out a trajectory on the unit sphere in 3-D generation space starting 
from the "high-energy" fixed point (1, 0, 0) at infinite scale and ending at the 
"low-energy" fixed point 4=(1, 1, 1) at zero scale. 

The fact that the mass matrix rotates with changing scales means that 
the usual definition of the 3 generation states as mass eigenstates has to be 
refined since it will have to be specified at what scales they are to be so 
defined. An analysis of the situation, defining each state at its own mass 
scale, leads to the following conclusion [TP] . For each fermion species T 
(i.e. whether U, D, L or N) the state vectors for the 3 generations form an 
orthonormal triad in generation space which is given in terms of the rotating 
vector r(/x) as follows: 

vi = ri, 

i"! A (ri A r 2 ) 



v 2 



|ri A (n Ar 2 )|' 



1-1 Ar " 2 (A\ 

where = r(mj) is the rotating vector r(/i) taken at the scale [i = mi, with 
rrii being the physical mass of the ith generation labelled from the heaviest 



1 There was an error in the coefficient of the right-hand side of eq. as given in ^Q, 
which means that the numerical value of p given in 0] and [H] should be multiplied by a 
factor of -y/5/3. Other results in these references are not affected. 
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(1) to the lightest (3). The elements of the mixing matrix (whether CKM 
! for quarks or MNS ^2] for leptons) are then given as the scalar products 
between the up-states v, and the down states v'-, thus: 




(5) 



Furthermore, for the fermion species U, D and L, the mass of the zth gener- 
ation state is given as the solution for to the equation: 



This criterion for determining the masses of the 3 generations, however, does 
not apply to neutrinos which, because of a likely see-saw mechanism [TB*] . 
may have physical masses different from the Dirac masses appearing in the 
above equation and hence require a special treatment to be explained later. 
With the above formulae, once given the trajectory for r, mixing matrix 
elements and mass ratios between generations can be evaluated, excepting 
for the moment where neutrinos are involved. 

The trajectory for r(p) depends on only 3 parameters, two of them corre- 
sponding to the vacuum expectation values of the Higgs fields which specify 
the rotation trajectory, and the third being the Yakawa coupling strength 
p in (J2J) which governs the rotation speed. The first 2 parameters are inde- 
pendent of the fermion species T, while the third p can in principle depend 
on T but, for consistency with the above prescription for defining masses of 
the generations, it was found that p has also to be T-independent (see jl] 
and later). In [S], these 3 parameters were fitted to m c /mt,mn/m T and the 
Cabibbo angle, which then allows one to calculate the whole trajectory via 
(J2J). The result is shown in Figure Q where from the location of the various 
states each marked at its own mass scale, the rotation speed with respect to 
p can be gauged. From this, the CKM matrix for quarks was calculated using 
(JSJ), and gave, as already mentioned, all elements within present experimental 
limits. 

2 v Oscillations 

Neutrinos require a special treatment because of the see-saw mechanism 
which is likely to give them physical masses different from the Dirac mass 
appearing in the mass matrix (0). On the one hand, the state vectors for 
neutrinos, in parallel to other fermions, ought to be defined at the scales of 



rrii = m T |rvVi 



(6) 
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their physical masses, while on the other the equation (JBJ) applies only to 
states whose physical and Dirac masses coincide, and hence not directly to 
neutrinos. It would appear therefore that some additional assumption will 
have to be made on the relationship between the physical and Dirac masses of 
neutrinos before the above scheme can be applied, and this was the approach 
we took in |14j . 

It turns out, however, that as far as reproducing the MNS mixing matrix 
elements is concerned, which is at present our main conern, then no such new 
assumption is needed. This arises as follows. The electron neutrino mass is 
now restricted by e.g. tritium decay experiments to below 3 eV jH], while 
oscillation experiments on solar an d atmospheric [7] neutrinos limit the 

mass differences between generations to less than 0.05 eV, implying that all 
3 generations have masses of at most order eV. According to the calculation 
cited in Figure ^ however, at scales of order eV, the vector r is already so 
close to the low energy fixed point ^(1, 1, 1) as to be indistinguishable in 
the figure. Referring next to the definition of the state vectors in one 
sees that in this situation when all 3 generations are close together, the 3 
vectors in the triad become respectively just the position vector, the tangent 
vector to the trajectory, and the vector normal to both the above, all taken 
at the same point. In other words, one has in particular for the state vector 
of the heaviest neutrino u s just: 

V! = v^=ro= -^=(1,1,1). (7) 

(Notice that neutrinos are conventionally labelled in the opposite order to 
that adopted in (jlj) above, i.e. from the lightest (z^) to the heaviest (^3)). 

Charged leptons being ordinary Dirac particles, their state vectors can 
be calculated by the method detailed in the preamble above in the same 
manner as for quarks. The result of the calculation performed in [S] and 
cited in Figure [T] gave: 

\r) = (0.996732,0.076223,0.026756), 

\n) = (-0.075925,0.774100,0.628494), 

\e) = (0.027068,-0.628482,0.777354). (8) 

From these and (J7J) above, using (j5J), one easily obtains the MNS mixing 
matrix elements: 

U l a = ( f M\v 3 ) = 0.7660, 

Uea = (e\v 3 ) = 0.1016, (9) 
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which are seen to be within the experimental limits obtained for these quanti- 
ties by respectively the atmospheric [7] and reactor |H] neutrino experiments: 

|C/ M3 | ~ 0.56 - 0.83, 

|C/ e3 | ~ 0.00-0.15. (10) 

One notices that the above result is independent of any assumptions about 
the Dirac masses of neutrinos, in contrast to the previous treatment of [T3] . 
In particular, because of the assumptions made then on neutrino masses, the 
treatment in ^3] is valid only for the range of masses permitted by the "vac- 
uum oscillation" solution to the solar neutrino problem. The present result 
in ©, on the other hand, having been derived independently of assumptions 
on neutrino masses, is valid for any solution to the solar neutrino problem, 
in particular for the large mixing angle MSW jTB] solution. This is of great 
practical significance, given the increasing preference of recent data for the 
LMA solution to the exclusion of the others [T71 [TH] . 

In principle, one can use the trajectory calculated in j3] and displayed in 
Figure ^ to evaluate also the tangent and normal vectors at the low energy 
fixed point r = ^(1, 1, 1), namely the state vectors and v n , and hence 
all the other elements of the MNS mixing matrix, in particular the solar 
neutrino angle U e 2- However, as we shall see in the next section, due to the 
limitations of the 1-loop calculation so far performed, such results would be 
unreliable. 

3 Limitations of 1-loop 

The 1-loop approximation with the parameters fitted as above and as in [3] is 
valid only near the fixed point (1, 0, 0) where the rotation is slow. Away from 
the fixed point, the RGE will receive large logaritmic terms from higher- loop 
contributions. As seen in Figure d however, down to the top mass scale at 
175 GeV, the vector r has rotated from the fixed point at infinite scale by 
an angle of only about 0.012 radians, while even from the top mass scale to 
that at the muon mass at 105 MeV, the vector has rotated by only about 
another 0.3 radians. For scales above the muon mass, therefore, one can hope 
a 1-loop approximation to have some validity. But at scales much below the 
/i mass, the 2-loop correction, which is expected to be of the order of the 
square of the 1-loop result, will become increasingly significant, thus making 
the 1-loop result unreliable. 
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Given then the limited validity range for the 1-loop calculation so far 
performed, one has to re-examine the results obtained to ascertain which 
could be regarded as reliable. Let us accept for present purposes the 1-loop 
approximation as sufficiently accurate for scales down to about the muon 
mass, in which range would then be included the 2 heavier generations of 
the 3 fermion species U, D, L, i.e. all species except the neutrinos. We 
notice from (J3J), however, that for U,D,L, the triad of state vectors for all 
3 generations are in fact already determined in this range, the vector v 3 for 
the lightest generation being given as the normal once the vectors for the 2 
heavier generations are specified. Hence, given the triads for both the U and 
D quarks, one can evaluate the whole CKM matrix with confidence. Indeed, 
this is born out by the comparison given in |3] of the calculation result: 



/ 


v ud 






\Vub\ \ 




v cd 




v cs 


\Vcb\ 


V 


V td 






\Vtb\ ) 



( 0.9745 - 0.9762 
0.217-0.224 



0.217-0.224 
0.9733 - 0.9756 
0.0486 



0.0043 - 0.0046 
0.0354- 0.0508 
0.9988 - 0.9994 



\ 0.0120- 0.0157 0.0336 
with the then available experimental numbers quoted from the databook [H] : 



/ 0.9745 
0.217- 
V 0.004 - 



- 0.9760 
0.224 
0.013 



0.217-0.224 
0.9737- 0.9753 
0.035 - 0.042 



0.0018- 0.0045 
0.036 - 0.042 
0.9991 - 0.9994 



(12) 



which are seen also to be within the most recent limits ^H]- On the other 
hand, the masses for the lightest generation, according to (JBJ), requires run- 
ning the vector r down to the mass scales of the lightest generation, and 
hence cannot be achieved reliably by the 1-loop approximation. And indeed, 
the value obtained in [3] for the electron mass was 6 MeV, to be compared 
with the correct empirical value of 0.5 MeV. The mass values for the light 
quarks u, d were also not well reproduced, but that may be due to intrinsic 
ambiguities in defining light quark masses. 

Next, turning to the MNS mixing matrix for leptons, we note first that 
the state vectors for all 3 charged leptons r, /x, e as quoted in (jSJ) ought, 
according to the above criterion, be reliable. The same, however, cannot be 
said for the state vectors of neutrinos, the mass scales of which are all way 
outside the range of validity of the 1-loop calculation. The only exception 
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is the vector for the heaviest neutrino v 3 which, as noted in the previous 
section, can be identified to a good approximation with the low energy fixed 
point r = 1, 1) known to be valid beyond 1-loop [20] • In consequence, 
one obtains the result © in good agreement with experiment. The other 
2 vectors in the neutrino triad both depend on the tangent vector to the 
trajectory and hence cannot be reproduced by the 1-loop calculation. Indeed, 
the value obtained in [3] for the solar neutrino angle U e 2, which depends on 
the state vector for the second heaviest neutrino u 2 , is about 0.23 and lies 
outside the experimental limits of about 0.4 to 0.7. 

In other words, the results of the 1-loop calculation on fermion mass 
and mixing parameters have now all been checked to agree with present 
experiment in the range they are expected to be valid, and are seen to deviate 
from data only outside the expected validity range. 

It should be emphasized, however, that this apparent success achieved by 
calculating only the diagram with a single (dual colour) Higgs loop makes 
sense only if one regards both the mass hierarchy and mixing as consequences 
of the mass matrix rotation. As far as the orientation of the rotating vector 
r(/x) rotation is concerned, it was shown in 0] that radiative corrections due 
to standard model particles are zero, and those due to other dual colour 
particles are small, leaving thus the calculated (dual colour) Higgs loop with 
no competition. The same approximation, however, would not be applicable, 
for example, to the normalization of the rotating vector. 

4 Analytic Approach 

Previous calculations with the DSM scheme, e.g. [5J, were done numerically, 
on which calculations the above remarks also rely. It is found, however, that 
with some simple yet quite reasonable approximations, an analytic approach 
is possible which provides first a check on the numerical results, and second 
an easier means for the reader to scrutinize and verify them. It also gives 
algebraic relations between fermion mixing elements and mass ratios which 
may be useful for certain purposes. 

First, let us parametrize the vector r(/i) in terms of the usual polar co- 
ordinates, thus: 

r(/i) = (cos 9(fi), sin 6(n) cos 0(/i), sin 6*(/i) sin 0(/x)), (13) 

and denote by <pi the corresponding angles taken at the scale /i = rrii with 
rrii being the mass of the state i. Next, we recall from the last section that the 
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1 1 

6»i (- cos(0 2 -4>i),-^- sin(0 2 - 0i) sin 2 , — sin(0 2 - 4>i) cos 2 ) 

t>2 &2 



1-loop approximation is valid only in the region near the high energy fixed 
point (1, 0, 0) so that there is no loss of generality at 1-loop to make a small 
angle approximation, i.e. assuming that all 0's are small and that, although 
the 0's are in principle arbitrary, their differences are small. In the same 
spirit, we assume also that 9\ is small compared to 02- In this approximation 
then, one has: 

Q2 

Vi = (1 - y, 0i cos 0i, 0i sin <M, 
v 2 = (0, cos 02, sin 02 ) 

+01 (- COs(0 2 - 

v 3 = (0, -sin 02, cos 2 ) 

+0isin(0 2 -0i)(l,-— cos02,-— sin0 2 ). (14) 

Using (fPfj) and (0), one obtains the approximate forms of the CKM matrix 
elements: 

e t )\ 

) + sin(0 s - C ) I j- sin(0 c - (f) t ) - j- sin(0 s - 6 ) 1 , 



cos(0 s - C ) I j- sin(0 s ~(f>b)-j- sin( 



!>s - 0fe) + 9 t cos(0 s - 0i), 

- 6 ) - 9 t cos(0 c - 0f), 

- 4> b ) - 9 t sin(0 s - (f> t ), 

» c - 0b) + #t sin(0 c - 0t). (15) 
From these formulae, one deduces immediately that: 



v* = 




v cs = 


cos(0 s - 


v ud = 




v cd = 


sin(0 s - 


v us = 


-v cd , 


V ts = 


—9 b cos( 


v cb = 


6 b cos(0, 


v td = 


6 sin(0< 


v ub = 


—0ft sin( 



\v*\ 


~ cos A0 


~1, 


\v cs \ 


= Kd ~ 


cos(A0) ~ 1, 


\Vcd\ 


= |Ks|~ 


— sin(A0) ~ small, 


Ws\ 


= 


A0 ~ small, 


\Vtd\ 


= 


0sin(A0) ~ very small 



(16) 
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In other words, simply on the premises of a rotating mass matrix as given 
in JTJ without using even the evolution equation (J2J) for the rotating vector 
r apart from the condition that r should remain near the high energy fixed 
point (1,0,0), one has already derived the well-known, and long-wondered 
at, hierarchy of CKM elements as observed in experiment. This result was 
anticipated already using some elementary differential geometry in [21] but 
is here now made completely explicit. 

Next, let us rewrite the evolution equation (j2J) for the unnormalized vector 
(x, y, z) as: 

-— = h ~— = h ~— = h (17) 

x dt ' y dt ' z dt ' 

with t = In//, k — 3p 2 /(327c 2 ), where 

fi = k ( + 4^4) . c y clic - ( 18 ) 

\x z + y A x z + z A J 

To get rid of the normalization, we take differences of these equations, so 
that we can still use the parametrization (|13p. which, with arbitrary and 
small 6, becomes: 

r ~ (1,0 cos 0,0 sin 0), (19) 

and 

fx ~ 2k, f 2 ^ /c(cos20 - 1), f 3 ~ -k(l + cos 20). (20) 
Hence we get for fx — fs and f 2 — f 3 respectively: 

1 d(8 sin 0) 



sin dt 
d{9 cos 0) 1 d(0sini 



cos dt sin dt 

which simplify to 

1 d(p 
— = —i 

sin 20 cos 20 dt 

Idd 



fc(3 + cos20), (21) 
2£;cos20, (22) 

(23) 



giving 



tan 20 
tan 20o 

6 
0o 



-— = -fc(2 + sin^20), (24) 



(25) 
(26) 




/ fj, \ ( cos 200 N 1 
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Equation (|26|) reduces to 

° ( ^ (27) 



6»o V n ) 

if we further assume that <fi 2 is small. 

With this solution, one easily checks that, for consistency with p has 
to be the same for all fermion species T to a good approximation, a result 
first found in [3] numerically, as follows. From (jUJ), one has: 

m 2 = m 1 sm 2 (8 2 - 6 1 ), (28) 

or else, recalling Q 2 ^> 6%, which is an approximation numerically similar to 
that for obtaining the linearized solution ()27|) above: 

m 2 ~ mi0f. (29) 

Hence, using fl2Tj) . one has: 

/^~»,(^f. (30) 



mi \m 2 



with 6>/ being the value of 9 at some chosen (large) initial value /i/. Applying 
the formula (J30|) to successively the fermion species U, D, L and eliminating 
6j, one obtains: 



f m c m b 
\mtm s 

( m c m T \ 
\m t m^ ) 



1/2 - x ^fcr, /„_ \ 2fc D 





Pi 




m c , 








Pi 











(31) 



which is consistent for arbitrary fij only when ku = = k^-, as required. We 
note that this result is independent of the choice of the 3 model parameters. 

Turning next to parameter-dependent results, we shall first determine the 
parameter values by fitting them to the 3 best known empirical quantities 
as was done in j^J, namely the 2 mass ratios m c /m t and m At /m T , and the 
Cabbibo angle. From (|3T!|) using the t scale as J, one deduces the approximate 
relations: 



m, 



■lk 



9 V ' 

m t \m c J 



Ak 



= ■ (32) 



ii 



Inputting then the empirical values in GeV of m t ~ 175, m c ~ 1.25, m T ~ 
1.777, m M ~ 0.105, one easily obtains; 



k ~ 0.21, ^~ 0.011. (33) 

These values compare very well with the values k ~ 0.20, 8 t ~ 0.012 obtained 
before from the numerical fit of jSJ , affording thus a check on both the present 
approximation and the previous numerical result. The good agreement is 
actually a little fortuitous, being beyond what can be expected from the 
crudeness of the approximations made in deriving the formulae used above. 

With now k and 6 t as 2 of our parameters, we can next proceed to evaluate 
other mass ratios and mixing matrix elements which depend only on 9. For 
mass ratios within the range of validity, we have only one relation (from (|31|) 
using nib = 4.2 GeV): 

x mfe y/(4 fe+ i) 

m s ~ m c [ — , (34) 

\m t J 

giving m s ~ 160 MeV (cf. 75 to 170 MeV [Hj). For the CKM matrix elements, 
we have: 

^l-^-l-i^)*, (35) 
giving V tb ~ 1 - 0.001 (cf. 0.9990 to 0.9993 0), and 

v u = -Va ~ -e b + e t ~ (^-) 2k , (36) 

y m t \rribJ 

giving V ts = —V c b ~ 0.05 (cf. 0.035 to 0.043 [6 ). These number also agree 
well with those obtained from the numerical result of [H] quoted in (|11|) . 

Inputting next the Cabbibo angle V us ~ 0.22 to determine the last pa- 
rameter which we can take as sin(0 s — <ft c ) and which for small <p differences, 
according to (fTB^l . can be take approximately as: 

sin(0 s - 4> c ) = — — (37) 
1 + (m s /m b ) 

we obtain sin(0 s — <f) c ) ~ 0.18. With this value fixed, we can now evaluate 
the remaining CKM matrix elements, thus: 

V ud ~ -V cs ~ 1 - l -V^ (38) 
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giving V ud Vcs ~ 0.976 (cf. 0.9734 to 0.9749 0), 



K 



V td ~ 6 sin(0 s - C ) ~ /_i f -£ j " (39) 

giving I4 d ~ .009 (cf. 0.004 to 0.014 6 ), and 



V ub ~6 b — sm{fa-fa)~ J — — — — — - ^, (40) 

t7s V rrit \m b J \ + (m s /m b ) 

giving ~ 0.004 (cf. 0.002 to 0.005 0). Again, these estimates agree 
well with the values (|TT| obtained before numerically in [Sj. Contrary to 
usual belief the CKM matrix elements appear related to the ratios of the two 
heaviest generations, which is expected if both come from a rotating mass 
matrix. 

Finally, using the values of the parameters determined, one can calculate 
also from ()14j) the 3 vectors v>, v M and v e for the charged leptons, as explained 
in the section above. We notice in ()14j) . however, that these vectors depend 
on the angles fa which may not be small, not just on their differences. These 
can be determined by solving equation (j23j) and using the result sin(0 s — fa) ~ 
0.18. This gives 2 solutions, cf) s ~ 0.61, 0.35, of which the fit in |3] correspond 
to the former. With the values 9 t ~ 0.011 obtained before and <p s ~ 0.61, 
one can then evaluate the corresponding angles for r and /i, and hence the 
vectors v T , v M , v e . One obtains: 

\t) = (0.994,0.070,0.029), 
\fj) = (-0.074,0.757,0.657), 

\e) = (0.0188,-0.657,0.757), (41) 

not far from the values quoted in JSJ) from jH] , and corresponding to the MNS 
lepton mixing matrix elements: 

U^ = (ji\vs) = 0.77, 

U e3 = (e\u 3 ) = 0.07, (42) 

which is seen to compare well with © and (jlOj) above. 

Hence, one sees that all the previous results obtained numerically before 
together with their agreement with experiment have now been confirmed by 
analytic considerations although in a rather crude approximation. 
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In summary, we conclude that a closer examination of the tenets of the 
DSM scheme has removed a previous restriction on its predictions on neutrino 
oscillations, making them now consistent with the favoured LMA solution for 
solar neutrinos. It has also explained away the couple of discrepancies noted 
before as limitations of the 1-loop calculation so far performed, which is 
seen otherwise to be in full agreement with experiment within its perceived 
range of validity. Furthermore, all previous numerical results have now been 
confirmed by analytic considerations. 
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